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Abstract 

The response of a CdW04 crystal scintillator to protons, a particles, Li, 
C, O and Ti ions with energies in the range 1-10 MeV was measured. 
The non-proportionality of CdW04 for low energy electrons (4 - 110 keV) 
was studied with the Compton Coincidence Technique. The energy depen- 
dence of the quenching factors for ions and the relative light yield for low 
energy electrons was calculated using a semi-empirical approach. Pulse-shape 
discrimination ability between 7 quanta, protons, a particles and ions was 
investigated. 
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1. Introduction 

Numerous astronomical and cosmological observations since 1930's sug- 
gest that most of the matter in the Universe is non-luminous and non- 
baryonic, while usual matter constitutes only ~ 4% of the Universe and 
the main components are dark matter 23%) and dark energy 73%) 
QQ. Dark matter (DM) is preferably related with particles which are neu- 
tral and only weakly interacting with matter (Weakly Interacting Massive 
Particles, WIMPs). Interaction of such particles with usual matter could 
be detected through observation of nuclear recoils created after scattering 
of WIMPs on atomic nuclei in sensitive detectors placed in low background 
conditions deep underground. Positive evidence of WIMPs observation at 
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~ 9a confidence level is reported in the DAMA experiment after 13 years 
measurements with large mass low background Nal(Tl) scintillators. This 
result was recently supported by the CoGeNT jl] and CRESST [5|] data. 
However, many other searches for WIMPs to-date gave only negative results 
[g], [^. Because of the fundamental importance of the question of the DM 
constituents, searches for WIMPs are under performance or planned in near 
twenty experiments, in particular, in the ton-scale projects EURECA with 
cryogenic Ge and different crystal scintillating bolometers [sj , DARWIN with 
scintillating noble gases ^ and DAMA/lton (proposed since 1996) [lOl] . 

When the energy released by the recoil ions is measured with scintillators, 
one of the main questions is the quenching of the scintillation light yield. 
Since a long time it is known that the amount of light produced in scintillating 
materials by highly ionizing particles is lower than that produced by electrons 
of the same energy yjj]. In a scintillator calibrated with 7 sources, signals 
from ions will be seen at lower energies than their real values, sometimes by 
more than one order of magnitude. Knowledge of quenching factor, QF (i.e. 
ratio of the measured ion energy in 7 scale to its real energy) is very important 
in searches for WIMPs and in predictions where the WIMPs signal should 
be expected. QFs depend on many factors such as: the scintillating material 
itself, its dopants and impurities, temperature; ion's Z and A numbers, ion's 
energy; time of collection of scintillating signal, etc. (see some examples in 
[1^). Since it is sometimes quite difficult to measure QFs for the needed 
ions and in the needed energy region (f.e. for low-energy heavy ions when 
the scintillation signal is expected at ~ 10 keV or less), some methods for 
the QFs estimation are of great interest. 

QFs for different ions in a scintillator could be not independent but related 
quantities. Such a h ypo thesis, supported by some experimental data, was 



discussed already in [11|. If true, on the basis of measurements of QFs for 
particles of one kind in some energy region (e.g. for a few MeV a particles 
from internal contamination of a detector), one would be able to calculate 
QFs for particles of another kind and for other energies (e.g. for low energy 
nuclear recoils). Further evidences in favor of this hypothesis were given in 
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Cadmium tungstate (CdW04) crystal scintillators are wi dely used in low 



counting experiments to search for 2/3 decay [13|, |l4, ll5|, |16|, |17[ and studies 



of rare (3 [18|, Il9| and a [20| decays. CdW04 has rather similar proper- 
ties to CaW04, ZnW04, CaMo04 and some other oxide crystal scintillators, 
promising targets for DM experiments 
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In this work we measured the response of a CdW04 crystal to few MeV 
energy ions and low energy (4 — 110 keV) electrons. Pulse-shape discrimina- 
tion ability was studied for protons, a particles, Li, C, O and Ti ions. Using 
the measured quenching for protons, we applied a semi-empirical approach 



12] to calculate QFs for other ions and low energy electrons. 



2. Crystal characteristics and signal processing 

The CdW04 crystal studied in this work is a parallelepiped of 10 x 20 x 25 
mm, manufactured by Scionix. CdW04 is a monoclinic, almost orthorhom- 
bic, crystal. From the specification of the manufacturer, the 10 x 20 mm 
faces are known to be cleavage (010) planes. Measurements performed with 
a diffractometei0 confirm this assignment, while the two faces 10 x 25 mm 
and 20 x 25 mm almost correspond to planes of indexes (100) and (001), 
respectively. 

The crystal is optically coupled through a 10 x 20 mm face to a pho- 
tomultiplier (PMT) ETL mod. 9256B. Its photocathode has an extended 
green response for a better matching with the light emission spectrum of the 
CdW04 crystal. In the measurements with low energy electrons the crystal 
has been wrapped with Teflon tape, while in measurements under beam an 
aperture of 10 x 10 mm has been opened in the tape covering one 20 x 25 
mm face, so that ions can enter directly. It has been verified with 7 sources 
that no measurable difference in light collection is found with or without 
aperture. The energy scale of the detector was established with ^^Na, ^''Co, 
^°Co, ^^'^Cs and ^^^Am calibration sources. 

In all the measurements the anode signals from the PMT are processed 
by a current to voltage converter which acts also as antialiasing filter and 
analysed by a transient digitiseiH already described in [21|] . The data, sequen- 
tially digitised every 50 ns, are stored in a temporary memory FIFO, which, 
in the presence of an event trigger, is stopped and read in a time interval 
starting about 30 /xs before the time of the event and extending to 128 jis on 
the whole. These data are stored on a mass memory for further analysis. In 
particular the amplitude spectra are obtained by the following procedure: for 
each event the mean value of the baseline is evaluated in the first 25 /is and 



^ These measurements have been performed at the "Centro di Cristallografia Strut- 
turale" of the University of Firenze. 

^The main characteristics are: 12 bits (ll-effective), 20 Ms/s, ±1024 mV Unear range. 
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subtracted channel by channel, afterwards the signal amplitude is obtained 
by summing up the channel contents in a time interval lasting 55 /is from the 
beginning of the signal. It is also checked whether the analysed waveform is 
in true or chance coincidence with the event trigger; overlapping waveforms 
are recognised and discarded. 



3. Response to protons, a particles, Li, C, O and Ti ions 

3.1. Measurements with pulsed beams 

Beams of protons, Li, C, O, and Ti ions, in the energy range 1 — 10 
MeV, were produced by the Tandetron accelerator of LABEC at the INFN- 



Florence. A dedicated beam line [22| provides short pulses by means of an 



electrostatic deflector which displaces the beam spot across a narrow window 
obtained by an adjustable slits system, housed at the entrance of the vacuum 
chamber, situated at the end of the beam-line and containing the crystal. 
The beam intensity and the width of the slits were adjusted to obtain the 
arrival of a single particle per pulse in a large majority of the cases. As 
already said in Section 2, in all these measurements the accelerated particles 
enter the crystal perpendicularly to a 20 mm x 25 mm (001) face. In these 
measurements the event trigger was derived from the voltage transition of the 
deflecting plates. Owing to this coincidence constraint on the stored pulse 
shapes, it is possible to build from them energy spectra almost free from 
background. The measured quenching factors for protons (1.0, 2.0, 3.0, 4.0 
MeV), Li (3.0, 4.5 MeV), C (3.0, 5.1, 7.5 MeV), O (6.0 MeV) and Ti (10.0 
MeV) ions are presented in Fig. 1. 

Differences of the pulse shapes, as a consequence of different specific en- 
ergy losses, have been also investigated. 

3.2. Measurements with a particles 

In addition to measurements with accelerated ions, the CdW04 crystal 
was irradiated with a particles. A collimated ^^^Am a source was placed, 
inside the crystal chamber, in front of the crystal (001) face and, at the same 
time, the crystal was irradiated by ^^'^Cs and ^"^00 7 sources in order to get 
a reference 7 energy scale. In this self-trigger was used and, taking 

advantage from the different zero-crossing times, a particle and 7 ray spectra 
were easily disentangled. From this measurement the quenching factor for 
the 5486 keV a particles, entering perpendicularly the (001) face, was found 
to be 0.127(1) as reported in Fig. 1. In a different series of measurements. 
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Figure 1: (Color online) Quenching factors for protons (•), a particles (■), Li (a), C (▼), 
O ("A") and Ti (0) ions measured with the CdW04 crystal scintillator. Solid lines represent 



calculations in accordance with 12|. On the right, lower QFs are shown in more detail. 



the three faces of the bare crystal have been irradiated with a particles. 
The measured QFs were: 0.197(1) for (010) face, 0.167(1) for (100) face and 
0.134(1) for (001) face. These measurements can be helpful for determining 
the response of the crystal to particles created in its volume and moving 
in random directions with respect to the principal planes. A comment is 
required by the small (< 5%) but significant difference between the two 
QF values relative to the (001) face: most probably, the difference could be 
attributed to the different light losses associated to the different locations of 
a and 7 light emitting sources (a surface point-like light source in the case 
of a particles and a volume distributed light source in the case of 7 rays) 
with respect to the different reflecting and diffusing properties of the bare 
and Teflon wrapped crystal. 

Difference in QFs for a particles moving in different directions, in addition 
to CdW04 detectors (here and e.g. in 20]) was observed also for ZnW04 [23[ 
and MgW04 2J] crystal scintillators. As quenching depends on ionization 



density [11|, such a behavior of the QFs has to be related with difference 
in ionization density and stopping power for a particles moving in different 
directions inside an anisotropic crystal. This interpretation is also supported 
by the behavior of the shape indicator (Section 3.4), which also depends on 
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the direction of the a particles, being closer to that of 7s when the quenching 
factor is larger. Differences of light collection due to surface effects could, in 
fact, change the quenching factor but hardly modify the shape of the light 
pulse. 

3.3. Calculations of quenching factors for ions 

Quenching factors Qi for protons, a particles and Li, C, O, Ti ions in the 



CdW04 scintillator were calculated following the work [12|, which is based 



on the classical Birks formula , as the ratio of light yield of an ion to that 



of an electron of the same energy: 

Qi{E) = Li{E)/L,{E), (1) 



where 




Here {dE/dr)i and (dE/dr)^ are the total stopping powers for ions cal- 
culated with the SRIM code |25| and electrons calculated with the ESTAR 
code |26[, respectively. 

Supposing the normalization factors Si^e equal for electrons and ions and 
independent on energy, we obtain that Qi{E) depends only on a single pa- 
rameter kB (the Birks factor). It is assumed that kB is also independent of 
energy and has the same value for all ions, if all data are measured under 



the same experimental conditions and treated in the same way [12[ . 

To determine the kB value for the present CdW04 measurements, the 
experimental points for protons were fitted by curve calculated with Eq. (1)- 
(3). The result is kB = 17.4 mg cm~^ MeV~^, and the obtained curve is 
shown in Fig. 1. Afterwards, QFs for all other ions were calculated with this 
kB value; all the results are also demonstrated in Fig. 1. 

As one can see, the calculations do not perfectly reproduce the experi- 
mental points, with the biggest deviation 30%) for Ti ions. Nevertheless, 
general agreement could be considered as acceptable, especially for a theory 
with only one parameter. Sometimes quenching factors in DM experiments 
are known with much bigger uncertainties, and calculations in accordance 
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with the above described approach could give valuable information on the 
expected QF values. 



3.4- Pulse-shape discrimination between 7 quanta, protons, a particles and 
light ions 

To test the ability of pulse-shape discrimination between 7 quanta {f3 
particles), protons, a particles and the light ions, the data were analyzed 
by using the optimal filter method proposed by E. Gatti and F. De Martini 
|27j (see also ^] where the analysis was developed for CdW04 crystal scin- 
tillators). For each experimental signal, its shape indicator was defined as 
SI = f{tk) X P(tk)/ Yl f{tk), where the sum is over time channels k, start- 
ing from the origin of signal and up to 60 /is, f(tk) is the digitized amplitude of 
a given signal. The weight function P{t) is defined as: P{t) = f^it) ~ f-yi^): 
where the reference pulse shapes fait) and f^(t) are the average of a few 
thousands shapes for a particles and 7 quanta, respectively, collected in the 
measurements with 7 and a sources. 

Distributions of the shape indicator versus energy for 7 quanta, protons, 
a particles, Li, C, O and Ti ions measured with CdW04 crystal scintillator 
are presented in Fig. 2. There is clear discrimination between 7s, protons, a 
particles and ions. The shape indicator values for protons and ions lie near 
the ones of a particles. This fact, however, does not necessarily imply that 
the pulse shape be equal for heavier ions and protons. 

It is worth to note also that the measurements with ^'^^Am a particles 
entering into different faces of the CdW04 crystal (Section 3.2) show a clear 
relation between quenching factor and shape indicator: for bigger QF^, when 
a signal is closer to signal caused by 7 quanta, SI^ value is lower, i.e. also 



closer to SI^. This effect was already observed in |20l . |23 



4. Measurements with low energy electrons 

It should be noted that the quenching of scintillation signals is observed 
not only for ions but also for low-energy (less than a few hundred keV) 7 
quanta and electrons. However, the nature of the electron-induced light- 
emitting excitation can have a different origin from that induced by heavy 
ions. Possible relation between QFs for ions and QFs for electrons is interest- 
ing and important as giving an additional method to obtain QFs for recoils 
induced by WIMPs in scintillators. 
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Figure 2: (Color online) Distributions of shape indicator (see text) versus energy for pro- 
tons, Li, C, O and Ti ions measured with the CdW04 crystal scintillator at the accelerator. 
Data for 7 quanta were obtained from a background run (mainly environmental radioac- 
tivity), while a particles were accumulated with the ^^^Am source installed in the vacuum 
chamber of the accelerator. 



A lot of work has been dedicated, also recently, to the study of the non- 
linearity of the response of scintillating crystals, coupled to different types 
of photo-detectors (photomultipliers in many cases) ^ to low-energy 7 quanta 
and electrons (see e.g. p 



30, iMl l32|, |33|, |34, l35|, |36[ and works quoted 



therein). We have studied the response of a CdW04 scintillation detector 
to low energy electrons by using electrons created in the CdW04 crystal by 
Compton scattered 7 quanta. 
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4.1. Experimental set-up 

To avoid surface effects for low energy electrons, the CdW04 crystal under 
study is used as a Compton spectrometeiEI. This technique requires that the 7 
rays, scattered by the crystal at an angle 9 with respect to the direction of the 
incoming coUimated monochromatic beam, be detected in coincidence with 
the signals generated in the crystal by the corresponding scattered electrons. 
However, when the energy of the scattered electrons is small, typically less 
than 20 keV, it becomes difficult to achieve reliable information both in 
energy and in timing, particularly for crystals, such as CdW04, whose main 
light emission is characterised by a long decay time (~ 14 /is) [s^: in that 
case a small anode signal consists of a sequence of almost randomly spaced 
single electron responses. 

To overcome this problem, we have used the production in opposite di- 
rections of the two 511 keV 7 rays from the singlet positronium following 
the decay of ^^Na. This variant of the Compton Coincidence Technique 
makes it possible to obtain a good timing of the event trigger through the 
coincidence between one of the 7 rays observed by a Nal(Tl) counter and the 
other one observed in a high purity germanium (HPGe) counter, after been 
scattered in the crystal under study, and therefore to study the response to 
electrons down to very low energy also for CdW04. Moreover the collima- 
tion of the 7 rays, impinging on the crystal, is performed via coincidences 
with Nal(Tl), so avoiding the border effects of a physical collimator. The 
front-end electronics for the Nal(Tl) detector consists of a passive integrator 
of the anode current signals followed by a JFET input linear buffer feeding 
a semi-gaussian shaping amplifier (Ortec mod. 572: time constant 0.5 /xs) 
whose bipolar output is sent to a Timing Single-Channel Analyzer (Ortec, 
mod. 551). HPGe signals, from one of the twin outputs of the charge pream- 
plifier, are processed in the same way as those of the Nal(Tl) detector and 
provide timing information, while the other preamplifier output is connected 
to a Gated Integrator Amplifier (Ortec: mod. 973) with integration time 5 
/is, and used for energy measurement. 

The geometry of the experimental set-up is shown in Fig. 3. The ^^Na 
source has nominal activity 400 kBq. ^^Na decays to the 1274.5 keV, first 
excited level of ^^Ne, through a [E^ax = 546 keV) annihilating in a 2 mm 
diameter zone. 



''This method was first proposed by Valentine and Rooney [371,1381]. 
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Figure 3: (Color online) Set-up for measurements of CdW04 crystal scintillator response 
to low energy electrons. All sizes are in mm. The angle 6 has been set to 12.5, 20 and 
35°. The Al cap surrounding the CdW04 crystal is not shown. 



As shown in Fig. 3, the source is on the axis of a system of two detectors: a 
3x3 Nal(Tl) (manufactured by Scionix) and the CdW04 crystal described 
in Section |2l The 10 mm thickness of this crystal causes an attenuation of 
the 511 keV 7 rays of ~60%. Special care has been devoted to minimise 
the material surrounding the CdW04 crystal and thus the probability of 
absorption and scattering along the path of the 511 keV 7 rays: to this 
purpose the cap surrounding the CdW04 is a cylinder of aluminium of only 
0.4 mm thickness. The 800 mm distance between source and front face of 
the Nal(Tl) detector defines a solid angle of maximum angular aperture of 
5.4°, so that a 511 keV 7 ray impinging on this detector has a companion 
511 keV 7 ray hitting the CdW04 crystal inside a circle of 16 mm diameter. 

The 511 keV 7 rays scattered at angle 6 from the crystal are detected by 
an HPGe detector (Ortec: mod. GMX 30 P: diameter 60 mm, height 60 mm, 
relative efficiency 30% and resolution (FWHM) 1.9 keV at 1.33 MeV). The 
HPGe detector is placed at ~ 200 mm from the centre of the crystal and, 
in successive measurements, at angles 9 of (12.5, 20.0, 35.0)° with respect 
to the axis defined by the Nal(Tl) detector and the CdW04 crystal. In this 
way, the overall angular range covered by the HPGe extends, with variable 
efficiency, from 0.5° to 44°, corresponding to energies of the scattered 7 rays 
from 510 keV to 395 keV and a complementary energy range of Compton 
electrons inside the CdW04 from 1 keV to 116 keV. 

It is worth to stress that with this set-up the problem of the time mark 
associated to even very low energy events in the crystal is overcome, because 
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the event trigger is obtained from the coincidence of 511 keV fuU-energy 
signals detected by the Nal(Tl) and signals detected by the HPGe in the 
energy range (390 - 511) keV. A drawback of this method concerns the 
intrinsic energy spread of the 511 keV 7 rays, which brings the FWHM of 
the 511 keV hne to 3.0 keV (the intrinsic HPGe resolution at this energy 
is ?s 1.7 keV). 

For monitoring gain and baseline of the HPGe, a ^^^Ba source is placed 
near the HPGe detector and shielded by a 25 mm lead shield from the CdW04 
crystal. For the same purpose, a source of ^"^^Am (emitting 60 keV photons) 
is placed near the CdW04. Lead shields have been also introduced which, 
to some extent, prevent 511 keV (and to a lesser extent 1274.5 keV) 7 rays 
from the ^^Na source, from directly reaching the HPGe detector. 

At the arrival of an event trigger, the linear signal from HPGe is peak 
analysed (by a 100 MHz, 8k ADC) and recorded, together with the infor- 
mation on CdW04 "signal shape", on a dedicated PC. As an example, the 
relevant part (470 keV - 520 keV) of the HPGe spectrum corresponding to 
7 rays scattered from the CdW04 crystal in a 4-days run (corresponding to 
about one tenth of the total statistics) at a mean angle 9 = 12.5°, is shown 
in Fig. 4a. The broad structure ("bump") in the 7 spectrum is due to 511 
keV 7s Compton scattered in the crystal, which are in prompt coincidence 
with the event trigger. The out-of-coincidence spectrum (Fig. 4a) clearly 
shows a 511 keV full energy peak, which is followed (because of multiple 
scattering) by a small tail extending from the full-energy peak down to the 
Compton edge. Obviously, a distribution of the same relative size is present 
at the low energy side of each channel of the bump and the combined effect of 
all channels produces the low-energy tail extending well below the minimum 
energy of the scattered 7s. This effect certainly contributes, with increasing 
importance, to the lower-energy part of the bump. Owing to this reason, 
only the upper part of the bump has been considered in the analysis. 

In Fig. 46, the corresponding digitizer pulse amplitude spectrum of elec- 
trons in prompt coincidence is reported. The nominal energy scale ("appar- 
ent electron energy" E^) is normalised to the full-energy peak of 511 keV, 
corresponding to a full absorption of the 7 inside the crystal. Due to the 
non-linear response, the energy scale does not reproduce the true energy of 
low-energy electrons. 

The small rise of the electron spectrum at the lowest energies is completely 
accounted for by the background distribution (see the out-of-coincidence 
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Figure 4: (Color online) a: Spectrum of the 511 keV 7 rays scattered by the CdW04 
crystal, acquired by HPGe in prompt coincidence with the corresponding 511 keV 7 rays, 
detected by the Nal(Tl) crystal. The two spectra are normalised to equal time intervals; 
the latter is then multiplied by 20. b: Distribution of pulse amplitudes from CdW04 con- 
verted to "apparent electron energy" in keV triggered by the prompt coincidence between 
HPGe and Nal(Tl); and out of coincidence normalised to equal time intervals. All the 
spectra refer to mean angle 8 = 12.5° (see Fig. 3). 



spectrum in Fig. 46). 

Within the statistical uncertainties, the energy of electrons associated to 
a 7 ray of energy is Ee = 511 keV - E^. The bump in the distribution of 
E^, shown in Fig. 4a, as well as those at ^ = 20° and 35°, have been divided 
in bins 2 or 4 keV wide, and for each bin the distribution of electron signal 
amplitudes from CdW04 has been obtained. Examples of the amplitude 
spectra of electrons are presented in Fig. 5 for mean energies 5.1 keV and 
19.1 keV (from 2-keV bins), and 92 keV (4-keV bin). Also in this figure, the 
electron pulse amplitude has been converted in keV to give the "apparent 
electron energy" E^. For each amplitude distribution, two values of the 
mean energy are given in the figure: the true mean energy of the scattered 
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electrons {E^ = 511 keV - (-E-y), coming from the mean value {E^) of the 
7 energy of the selected bin on the 7 bump, and the mean of the apparent 
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Figure 5: (Color online) Energy spectra of electrons with energies 5.1, 19.1 and 92.0 
keV measured by CdW04 crystal scintillator. The solid lines represent the fit to the 
distributions by a gaussian function. 



Distributions, as those shown in Fig. 5, have been used for two purposes: 

1) from each electron distribution the value of FWHM has been extracted 
and compared with those obtained by irradiating the CdW04 crystal with 7 
rays of similar energies from radioactive sources; 

2) for each bin of the selected regions of the three 7 bumps (at 6 = 12.5°, 20°, 
35°) the relative light yield (Ef,/E^ of the corresponding electron energy 
distribution has been evaluated. 

4-2. Response of CdWO^ crystal scintillator to low energy electrons 

The dependence of the energy resolution on energy of electrons and 7s 
from radioactive sources in CdW04 is reported as a function of energy in Fig. 
6. The main findings are that no appreciable differences between electrons 
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and 7 rays are apparent in the overlapping energy region. The continuous hne 
shows a minimum-x fit of the data with the function FWHM(%)= ./oJE; 
(a = 38200±400 keV, x^/n.d.f. = 0.45, where n.d.f. is the number of degrees 
of freedom). 
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Figure 6: (Color online) FWHM of the distributions of pulse amplitudes from the CdW04 
counter for Compton electrons and 7s from various radioactive sources. 



Concerning the relative light yield, examples of distributions of the ratios 
Ee/ Ee, for the same energy bins of Fig. 5, are given in Fig. 7. From this type 
of distributions, and for each 4-keV bin in the HPGe bumps at angles 6 of 
(12.5, 20.0, 35.0)°, the relative light yield is derived as R = (eJE^Y They 
are presented as a function of the true electron energy in Fig. 8. Error bars 
include the statistical uncertainties as well as estimated systematic errors 
related to the measurement itself and to the procedure of data analysis. The 
relative light yield remains substantially stable at the 100% value above 80 
keV, decreases slowly at lower energies down to 95% at 20 keV and more 
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deeply below this energy. 
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Figure 7: (Color online) Examples of distributions of ratio E^/Ee of the CdW04 scintil- 
lation response together with their fits by a gaussian function. 



In CdW04, no indication is found of an increase of the luminescence ef- 
ficiency with decreasing electron energy, as reported for Nal(Tl). To our 
knowledge, no previous measurements with internal electrons have been re- 
ported for CdW04. Compared with Fig. 10 of [30], our Fig. 8 shows a slower 
decrease of the R for energies below 80 keV. We must remark, however, that 



the data of [30[ are obtained by excitation with 7 rays, and show in fact the 
well known anomalies at the K edge. 

We also calculate here the relative light yield as 

R = L,{E)/E (4) 

with Le obtained with Eq. (3) and the Birks factor kB = 17. 4 mg cm~^ MeV~^. 
the same as was used to describe the experimental data for ions (Fig. 1); R 
is normalised to 1 at 100 keV. However, as one can see from Fig. 8, the ex- 
perimental data for electrons are better described with the value kB = 9.0 
mg cm~^ MeV~^, when the points at > 20 keV are perfectly described. 
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Figure 8: Relative light yield (normalised to 1 at 100 keV) as a function of the electron 
energy. Lines show calculations with Eq. (4) with different kB values. 



while disagreement still is present at lower energies. One could expect that 
introduction of kB dependent on energy (see |4o[. l4l[ 42|) or introduction of 



additional terms in the Birks equation (see e.g. 43|) would allow to improve 
the description of the data. 



5. Discussion and conclusions 

(1) It should be stressed, the non-proportionality in the scintillation re- 
sponse is an important characteristic to estimate the energy resolution achiev- 
able with scintillation crystals ISO]. In accordance with current understand- 
ing (see recent reviews 33|, |3J, [SSi]), non-proportionality of the relative light 
yield for electrons and 7 quanta at low energies is one of the main rea- 
sons of poor energy resolution of a scintillator for 7s even at higher energies 
because of a high probability for a 7 quantum to interact with a detector 
more than once and to lose the total energy through creation of few lower 
energy electrons. Many scintillating materials were studied till now (see 



331 . 134 I35I . l36l . I44J . I45I . kg and refs. therein), and the best energy resolution 
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(better that 3%) was reached with scintillators which show good proportion- 
ality of the light yield down to ~ 10 keV (f.e. LaCl3(Ce), LaBr3(Ce), Srl2(Eu) 
fii']). While response of CdW04 to low energy 7 quanta was already studied 
[30. . |45|, results for low energy electrons (down to 5 keV) are presented here 
for the first time. 

(2) Measurements with the accelerator for protons, Li, C, O and Ti ions 
with energies in 1 - 10 MeV interval, together with the 5.5 MeV a particles 
from ^^^Am, allow to obtain data on quenching factors interesting for the 
dark matter studies. Description of these data with Eq. (l)-(3) with the 
same value of the Birks factor kB = 11 A mg cm" -2 MeV"^ (obtained by 
fitting the data for protons), while being not perfect, nevertheless shows 
relevant agreement, with the biggest deviation of ~ 30% for Ti ions. This 
further supports the hypothesis [12|| that quenching factors for all ions are 
not independent and could be described with the same kB value, if the data 
are collected in the same experimental conditions and are treated in the same 
way. This approach allows to obtain QFs for low energy nuclear recoils on the 
basis f.e. of QFs values for a few MeV a particles from internal contamination 
of a detector. Sometimes the QFs for low energy nuclear recoils are known 
with uncertainties much higher than 30%, and the description with Eq. (1)- 
(3) could provide an important estimation of the needed QFs. 

(3) It is possible to describe the obtained data on the relative light yield 
for low energy electrons in the CdW04 by Eq. (4) in a perfect way for E' > 20 
keV, however still with disagreement at lower energies. Similar description 
on the basis of the Birks equation was obtained previously for CaW04 in 47 



and for liquid scintillators in the Double Chooz experiment in |48| . However, 



for CdW04 the value kB = 9.0 mg cm~^ MeV~^ determined for such a 
description is different from that obtained from fit of the data for the ions: 
kB = 17 A mg cm~^ MeV~^. This unfortunately closes the additional way 
to obtain QFs for low energy nuclear recoils relevant for the DM searches 
from the light yield non-proportionality for the electrons measured with the 
same scintillator. This conclusion is not unexpected; it is supported also 
by the following considerations. As it is known (see f.e. |35|), for some 
scintillating materials (Nal(Tl), CsI(Tl)) the relative light yield for 7 quanta 
and electrons is > 1 at low energies. At the same time. Eg. (4) can describe 
only quenched values R < 1 (as those in Fig. 8 or in 47, Hsj) but not the 
enhanced R > 1 (giving at most only R = 1 with kB = 0), thus it is not 
suitable for Nal(Tl). One could assume that the mechanisms of quenching 
for ions and non-linear response to low energy electrons (7s) are different. 
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Nevertheless, description for ions by Eq. (l)-(3) is valid also for Nal(Tl); see 
f.e. Fig. 13a in 12| where the kB value obtained by fitting data for Na recoils 
allowed to perfectly describe QFs for I recoils. 
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